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摘 要 : 植被 覆盖 度 和 生产 力 指标 作为 评价 草原 生态 系统 退化 已 有 较 多 研究 ,各 单一 评价 指标 互相 独立 难以 对 不 


同 草 原 退 化 程度 进行 综合 评价 。 本 研究 以 新 疆 天 山中 部 巴 音 布鲁克 草 1 提出 了 一 种 基于 标准 化 处 理 


分 指标 耦合 的 草原 退化 遥感 评估 方法 。 选 取 草 原 植被 覆盖 度 . 草 层 平均 高 度 


总 产 草 量 3 个 分 指标 ,通过 主 成 分 分 


析 法 确定 分 指标 权重 ,并 引入 Min-Max 标 准 化 方法 处 理 并 构 建 综合 评价 指标 间 地 退化 指数 (Grassland Degradation 
Index, GDI) ,最 终 通 过 Landsat 影 像 反 演 以 及 草地 退化 指数 变化 率 的 合理 分 级 评估 了 1986 一 2021 年 新 疆 天 山 巴 音 


布鲁克 草原 的 退化 程度 。 结 果 表明 :(1) GDI, 与 NDVI 的 相关 性 最 好 。 


总 面积 的 比例 为 60.51 %; 不 同 草地 类 


(2) 2021 年 巴 音 布鲁克 草原 未 退化 面积 


型 的 草地 退化 程度 有 明显 差异 ;空间 分 布 上 表现 为 由 盆地 向 山地 退化 加 重 的 


趋势 ,(3) 通过 辐射 配 准 方 法 可 以 将 CDI 模型 应 用 于 其 他 年 份 ,2000 一 2009 年 巴 音 布鲁克 草原 退化 程度 明显 改 


善 ,2009 一 2021 年 草原 退化 程度 轻微 波动 。 研 究 结果 将 为 指导 巴 音 布鲁克 草原 退化 程度 的 评价 和 保护 草原 生态 系 


统 提供 一 定 的 数据 支撑 和 理论 基础 。 


关键 词 : 草原 退化 ; Landsat 影 像 ;， Min-Max 标 准 化 ; 草地 退化 指数 ; 


草原 是 我 国 陆 地 生态 系统 的 重要 组 成 部 分 , 占 
国土 面积 的 41.67% ,其 中 各 类 天 然 草 地 约 有 4x 
10° hm2 ,可 利用 面积 为 3.9x10 hm ,集中 分 布 在 中 西 
部 地 区 ”。 新 疆 是 我 国 五 大 草原 畜牧 业 省 份 之 一 ， 
天 然 草地 毛 面积 5.6x10 hm2, 可 利用 面积 4.8x107 
hm ,位 居 全 国 第 三 位 ,草地 面积 约 占 新 疆 总 面积 的 
34.44%5 。 然 而 , 受 自 然 和 人 为 因素 影响 ,我 国 草原 
正 发 生 不 同 程度 的 退化 "。 新 疆 天 然 草 地 已 有 85% 
出 现 不 同 程度 的 退化 ,与 20 世 纪 60 年 代 相 比 , 巴 
音 布 鲁 克 草原 草地 植被 覆盖 度 由 原来 89.4% 下 降 到 
30%~50%5 , 巴 音 布鲁克 草原 退化 面积 达到 30.48x 
10° hm , 占 草地 总 面积 的 19.61%53。2000 年 以 来 ， 
我 国 实施 了 包括 退 牧 还 草 和 草原 生态 奖 补 机 制 在 
内 的 一 系列 重大 生态 工程 ,重点 进行 退化 草原 修 
复 治 理工 作 。 


回归 模型 ;， 巴 音 布 鲁 克 草原 


技术 被 广泛 应 用 到 草地 监测 方面 。 草 原生 态 系统 
研究 中 ,学 者 们 大 多 选用 植被 指数 .植被 覆盖 度 、 植 
被 生产 力 等 指标 来 评价 草地 状况 。 吴 晓 全 等 ” 利 
用 CASA 模型 估算 了 新 疆 天 山地 区 的 植被 净 初 级 生 
产 力 的 时 空 分 布 并 分 析 其 驱动 因素 ,结果 表明 天 山 
地 区 植被 - 净 初 级 生产 力 空间 上 呈 西 高 东 低 .由 北 
向 南 递减 的 特征 ;针对 草原 退化 或 健康 状况 评价 问 
题 , 杜 自强 等 "选取 植被 覆盖 度 、 牧 革 地 上 生物 量 
和 牧草 可 食 率 3 个 指标 加 权 平 均 构建 了 草地 退化 指 
数 , 通 过 TM 影像 遥感 监测 黑河 中 上 游 草 地 1986 一 
2003 年 的 退化 情况 ,结果 表明 草地 退化 加 剧 、 局 部 
改善 .整体 恶化 的 变化 格局 ; 赵 玉 婷 等 基于 “ 压 
力 -状态 -响应 ”(PSR ) 模 型 ,从 压力 状态、 响应 3 个 
层面 选取 17 个 指标 构建 了 甘南 藏族 自治 州 高 寒 牧 
区 草地 生态 系统 健康 评价 体系 ,并 表明 13 a 来 该 地 


草原 幅员 辽阔 , 仅 采 用 野外 调查 手段 不 能 满足 
大 尺度 评价 要 求 , 且 监测 周期 长 工作 量 大 ,因此 3S 


收 稿 日 期 : 2022-06-06; 
基金 项 目 : 第 三 


修订 日 期 : 2023-01-28 


区 高 寒 草 地 健康 状况 负 向 变化 趋势 明显 ; 陆 均 等 
通过 CVOR 指数 研究 表明 , 巴 音 布鲁克 高 寒 草原 生 


次 新 疆 科 学 考察 项 目 (2022xjkk0401) ; 王 宽 诚 教育 基金 会 (中 亚 成 海 生态 与 环境 修复 国际 研究 团队 ) 


作者 简介 : 赵 剑 (1997-) , 男 ,硕士 研究 生 ,主要 从 事 草原 生态 遥感 监 疯 


通讯 作者 : ELE. E-mail: likh@ms.xjb.ac.cn 


636 — 646 页 


I£ ZS 
|. E-mail: 715171198@qq.com 


http: //azr.xjegi.com 


202305.00262v1 


chinaXiv 


态 系 统 健康 状况 由 2004 年 的 一 般 病态 逐渐 恢复 为 
2012 年 健康 状态 , 且 退 化 状况 由 强度 退化 恢复 为 中 
度 退 化 ; 吴 志 丰 等 ”探究 了 一 种 基于 遥感 的 参照 覆 
盖 度 提取 及 草地 退化 评价 研究 的 方法 ,在 一 定 程 度 
上 解决 了 目前 大 范围 草地 退化 评估 过 程 中 面临 的 缺 
乏 参 照 系 统 及 遥感 数据 应 用 误区 等 问题 ;李霞 等 ” 
选取 盖 度 .总 产量 和 可 食 草 产量 3 个 指标 并 结合 GIS 
中 的 空间 插值 法 反 演 甘肃 省 草地 退化 状况 ,结果 显 
示 甘 肃 省 退化 草原 面积 1.79x10 hm2?, 占 全 省 草原 总 
面积 的 69.65%。 草 原 退 化 程度 反映 在 植被 覆盖 度 、 
地 上 生物 量 .可 食 草 比例 草地 退化 指示 物种 以 及 
土壤 有 机 质 含量 等 指标 的 变化 上 ,诸多 研究 对 不 同 
指标 进行 了 分 析 , 但 退化 表征 指标 复杂 且 受 限于 各 
种 调查 技术 ,各 个 分 指标 相互 独立 ,综合 评估 指标 
的 建立 是 面临 的 问题 。 

目前 ,对 草原 退化 程度 的 分 级 方法 种 类 繁多 ， 
各 退化 评价 指标 差异 大 且 难 统一 ,缺乏 综合 评价 指 
标 体系 ,通过 计算 分 指标 的 权重 值 ,实现 不 同 指标 
的 有 机 综合 是 提高 草原 监测 精度 的 关键 ””。 因 此 ， 
本 研究 选取 植被 覆盖 度 .总 产 草 量 以 及 草 层 平 均 高 
度 3 个 指标 ,并 引入 了 Min-Max 标 准 化 方法 对 分 指 
标 进 行 处 理 ,进一步 构建 草地 退化 指数 (Grassland 
Degradation Index, GDI) ,以 1986 年 野外 调查 数据 作 
为 未 退化 参照 ,通过 中 高 分 辨 率 的 Landsat 影 像 反 演 
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并 计算 得 出 GDI 的 变化 率 来 表征 草原 退化 情况 , 旨 
在 建立 草原 退化 遥感 综合 评估 指标 ,解决 草原 退化 
大 尺度 监测 效率 低 的 问题 ,评估 20 世纪 80 年 代 以 
来 巴 音 布鲁克 草原 的 退化 情况 。 


1 材料 与 方法 


1.1 研究 区 概况 

巴 音 布鲁克 草原 位 于 新 疆 维吾尔 族 自治 区 巴 
Ef FB a h A Yea A pa EL a AL LB, AR PE 
292 km ,南北 宽 108 km, 主要 由 大 尤 尔 都 斯 盆地 和 
小 尤 尔 都 斯 盆地 以 及 天 山 山 脉 中 段 组 成 ,海拔 为 
1624~4606 m。 其 气候 为 典型 高 寒气 候 类 型 , 雨 雪 
天 气 多 ,冬季 漫长 夏季 短暂 ,年 均 气 温 -4.8 “C ,全 年 
积 雪 日 约 137 d, 积 雪 深 度 45 cm'”。 草 地 类 型 以 沼 
泽 类 .高寒 草 甸 类 和 高 寒 草 原 类 为 主 , 地 表 优 势 植 
WARS (Kobresia capillifolia ) RF (Polygo- 
num viviparum) Z ENF (Stipa purpured)、 寒 生 羊 
F (Festuca kryloviana ) 等 。 

地 面 调查 样 地 点 数据 共 52 个 ,其 中 建 模 集 样 点 
35 个 ,验证 集 样 点 17 个 (图 1)。 
1.2 数据 来 源 与 处 理 

本 文采 用 与 研究 相关 的 遥感 数据 和 地 面 调 查 
数据 。 其 中 20 世 纪 80 年 代 草 原 普查 数据 和 2021 年 
地 面 样 地 调查 数据 均 来 自 于 中 国 科 学 院 巴 音 布 鲁 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2019)3333 号 ,对 底 图 边界 无 修改 。 


图 1 研究 区 位 置 及 样 地 点 位 


Fig. 1 Location and sample points of the study area 
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克 草 原生 态 系统 研究 站 ;选取 的 遥感 影像 为 2000 年 
7 H 6 H Landsat ETM+、2009 年 7 月 31 H Landsat 
ETM+ 和 2021 年 7 月 24 日 Landsat OLI 共 3 期 影像 ， 
云 量 低 于 20% 且 与 实地 采样 时 间 最 为 接近 ,数据 间 
隅 控制 在 10 a 左右 ,保持 一 定 的 逻辑 ;数字 高 程 
DEM 数据 来 自 于 地 理 空间 数据 云 (https://www.gs- 
cloud.cn)。 

预 处 理 ; 

(1) 遥感 影像 (ENVI 5.3 软 件 ):Radiometric Cal- 
ibration .QUAC „Seamless Mosaic .边界 裁剪 。 

(2) 地 面 调查 数据 (ArcGIS 10.6) : Add Data T. 
具 添 加 坐标 点 ,Extraction Values to Points 模型 构建 。 

(3) 数字 高 程 数据 (ArecGIS 10.6) :拼接 ,边界 裁 
剪 和 格式 转换 。 
1.3 统计 学 方法 
1.3.1 回归 分 析 回归 分 析 (Regression) 是 一 种 试 
图 以 一 个 或 多 个 自 变 量 来 解释 另 一 个 因 变 量 的 统 
计 方 法 ,在 非 线性 回归 模型 中 ,一 种 类 型 是 可 以 通 
过 变量 变换 ,化 为 线性 模型 进行 分 析 , 即 曲线 回归 
(曲线 拟 合 ) ,一 般 形式 为 : 

y;=f (x,0)+e, (1) 

式 中 : y, 为 因 变 量 ( 被 解释 变量 ); f(x,0) 为 参数 佑 
计 方 程 ; e; 为 误差 的 函数 , 拟 合 后 取 最 小 值 。 
1.3.2 Min-Max 标 准 化 ”由 于 不 同 的 数据 有 不 同 的 
量 纲 ,而 进行 综合 指标 构建 时 需要 有 统一 的 标准 ， 
因此 需要 对 数据 首先 进行 标准 化 处 理 。 数 据 标准 
化 可 以 消除 不 同 数据 之 间 不 同 量 纲 的 影响 ,解决 数 
据 不 匹配 的 问题 ,常用 的 标准 化 方法 有 Min-Max 标 
准 化 (Min-Max Normalization ) 7 z-score 标准 化 方法 
等 ,考虑 数据 量 和 数据 处 理 的 便捷 性 ,本 研究 采用 
的 方法 为 Min-Max 标 准 化 方法 ,标准 化 处 理 之 后 数 
据 在 [-1,1] 范 围 内 。 


X, = 


X-X; 
Tom Xan i 
式 中 : 成为 归 一 化 处 理 后 的 数据 ; X 为 原始 数据 ; 
Xi 为 数据 的 最 小 值 ; 和 ,为 数据 的 最 大 值 。 

1.3.3 主 成 分 分 析 法 “指标 权重 常用 的 确定 方法 有 
BRAT ANE 、 主 成 分 分 析 法 “、 聚 类 分 析 法 ” 
等 。 但 专家 打分 法 摊 杂 了 太 多 的 主观 意识 , 主 成 分 
分 析 (Principal Component Analysis ,PCA ) ,是 一 种 通 
过 正 交 变换 将 一 组 存在 相关 性 的 变量 转换 为 一 组 
线性 不 相关 变量 的 数学 统计 方法 ,将 数据 的 大 部 分 


信息 以 贡献 程度 靠 前 的 变量 集中 表示 。 已 有 多 项 
学 者 的 研究 表明 , 主 成 分 分 析 法 应 用 于 草地 退化 分 
旨 标 的 权重 确定 具有 一 定 的 科学 性 和 准确 性 ””。 
1.3.4 精度 评价 ”GDI 模型 的 精度 采用 验证 点 位 的 
遥感 拟 合 值 与 地 面 点 位 实际 计算 值 的 均 方 根 误差 
结果 来 进行 衡量 , 均 方 根 误差 ( RMSE ) 是 用 来 衡量 
观测 值 同 真 值 之 间 的 偏差 , RMSE 对 相差 较 大 的 数 
据 相关 性 更 高 ,所 以 其 针对 异常 值 更 为 敏感 ,计算 
公式 为 : 


-lx 2 
RMSE = ES, -co1,) 


i= 


(3) 


式 中 : RMSE 表示 均 方 根 误差 ; n 表示 样本 的 个 数 ; 
GDI, 表示 地 面 调 查 点 位 实际 计算 值 ; GDI, 表示 遥 
感 拟 合 值 。 在 用 于 评价 反 演 模型 的 精度 时 , RMSE 
值 越 小 表示 模型 的 精度 越 好 。 

1.4 草原 退化 分 级 方法 

1.4.1 植被 指数 ” 归 一 化 植被 指数 (Normalized Dif- 
ference Vegetation Index,NDVI) 是 应 用 最 广泛 的 植 
被 指数 ,可 以 有 效 反映 植被 状况 , 且 计 算 简 单 ,研究 
也 最 深入 ”。 此 外 ,本 文 义 引入 了 其 他 植被 指数 ， 
差 值 植被 指数 (Difference Vegetation Index, DVI) 、 比 
值 植被 指数 (Ratio Vegetation Index, RVI) .增强 型 植 
被 指数 (Enhanced vegetation index, EVI) .土壤 植被 
指数 (Soil-Adjusted Vegetation Index, SAVI)( 表 1)， 
并 给 出 其 计算 公式 。 

1.4.2 草地 退化 指数 (GDI) 草地 退化 指数 (Grass- 
land Degradation Index,GDI) 是 大 尺度 上 用 于 评价 草 


表 1 植被 指数 及 计算 公式 


Tab. 1 Vegetation index and calculation formula 


植被 指数 计算 公式 参考 文献 
归 一 化 植被 指数 NDVI = Rx Rea 23] 
Ry, + Raa 
差 值 植被 指数 DVI=Ry, ~ Rio 24] 
HEATH RVI- As 25] 
增强 型 植被 指数 = Rye 7 Rr 26] 
BV I=22 X R OUR -7.5R +1 

土壤 植被 指数 SAVI= A A +L) 27] 
三 波段 梯度 差 植 ay Ry, = Ria _ Ri mi Re 28] 
被 指数 ete A Nir A Red A Red 一 À creen 


4 TGDVI <0, 则 TGDVI=0 


TE: Ra 表示 近 红 外 波段 ; Ria 表示 红色 波段 ; Ru 表示 蓝 色 波 段 ; 


RR。 表示 绿色 波段 ; Aw 表示 近 红外 波段 的 波长 ; A 表示 红外 波段 


的 波长 ; A 表示 绿色 波段 的 波长 ; L IERI RZ, L=0.5 。 
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植被 覆盖 度 .总 产 草 量 和 草 层 平均 高 度 3 个 因子 的 
加 权 综 合 , 是 草原 植被 状况 的 数字 化 表示 , 取 值 范 
围 为 0~1, 计 算 公式 为 : 

GDI, = Dow, (4) 


式 中 : GDI, 表示 地 面 调 查 点 位 草地 退化 指数 ; n 表 
示 指 标的 个 数 ; v, 表示 i 指标 的 标准 化 值 ; w 表示 i 
标的 分 权重 。 

143 草原 退化 分 级 标准 ”以 1986 年 草原 样 方 调 
查 资料 作为 未 退化 参照 标准 ,根据 《天 然 草 地 退化 、 
沙化 . 盐 渍 化 的 分 级 指标 》(GB 19377-2003) AIRT 
人 的 研究 2 ,参考 了 植被 覆盖 度 减 少 率 (未 退化 < 
10% 、 轻 度 退 化 11%~20% .中 度 退化 21%~30% . E 
度 退 化 >30% ), 草 层 平 均 高 度 减少 率 ( 未 退化 <10%、 
轻 度 退化 11%~20% .中 度 退 化 21%~50% 、 重 度 退 
化 >50% ) 和 总 产 草 量 减少 率 (未 退化 <10% 、 轻 度 退 
化 11%~20% .中 度 退 化 21%~50% 、 重 度 退 化 >50%) 
的 单 指标 分 级 阔 值 ,加 权 平 均 得 到 了 综合 指标 草地 
退化 指数 的 分 级 阔 值 ,将 研究 区 草原 退化 程度 分 为 
4 个 等 级 :未 退化 (GDI 的 减少 率 <10% ) 、 轻 度 退 化 
(GDI 的 减少 率 为 10%~20%) .中 度 退 化 (GDI 的 减少 
率 为 20% ~40.8% ) .重度 退化 (GDI 的 减少 率 > 
40.8% )。 


2 结果 与 分 析 


2.1 20 世 纪 80 年 代 草地 退化 指数 的 计算 结果 

草原 退化 的 参考 基准 值 根据 20 世 纪 80 年 代 草 
原 普查 数据 确定 ,新 疆 天 山 巴 音 布鲁克 草原 分 为 高 
KEAK ,高 寒 草原 类 沼泽 草地 类 3 种 草地 类 ,地 
表 植 被 包含 19 种 植物 ( 表 2)。 

结合 遥感 影像 和 野外 实测 样 地 精细 划分 29 个 
草地 类 型 ,计算 20 世纪 80 年 代 草 原 普 查 数据 形成 
了 退化 背景 的 GDI, 参考 值 ( 表 3)。 GDI, 最 大 值 出 
现在 高 宕 草 甸 类 的 线 叶 渍 草 、 杂 类 草 类 型 ,表明 该 
类 型 植被 状况 最 好 , GDI, 最 小 值 出 现在 高 寒 草原 类 
寒 生 羊 茅 . 冰 草 ,该 类 型 植被 状态 最 差 ,整体 来 看 ， 
高 寒 草 甸 类 和 沼泽 类 植被 要 好 于 高 寒 草 原 类 。 
2.2 2021 年 草地 退化 指数 的 计算 结果 

2021 年 地 面 调查 点 位 的 70% 作 为 建 模 集 样本 , 
共 35 个 样 地 点 位 , 即 高 寒 草 甸 类 样本 点 6 个 ,高寒 
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表 2 主要 植物 种 类 
Tab. 2 Main plant species 


草地 类 植物 种 类 
高 寒 草 甸 类 AREH (Carex oxyleuca ) 


鬼 箭 锦 鸡 儿 (Caragana jubata ) 
2K litte AL (Kobresia capillifolia) 
细 果 蔓草 (Carex stenocarpa ) 
ERE (Polygonum viviparum ) 
Asta A (Carex turkestanica) 
H (Koeleria macrantha ) 


SEE FF (Festuca kryloviana ) 


宽 穗 赖 草 (Leymus ovatus ) 

冰 草 (Agropyron cristatum ) 

RAW BSE Uris loczyi Kanitz) 
紫花 针 茅 (Stipa purpurea ) 

HUES (Artemisia dracunculus ) 
EZ (Secale sylvestre) 

19 (Artemisia frigida ) 
EAE ET (Stipa subsessiliflora ) 
WIEZE (Carex stenophylla) 
布 顿 大 麦草 (Hordeum bogdanii) 


Wi KT 2S EE (Juncus compressus ) 


沼 译 草地 类 


草原 类 样本 点 24 个 ,沼泽 类 样本 点 5 个 ( 表 4) ,采样 
时 间 为 7 月 10 一 31 日 ,植被 覆盖 度 主 要 通过 目 视 佑 
测 法 得 到 ,每 个 样 地 包含 3 个 1 mxl nm 的 小 样 方 , 通 
过 间 制 法 获取 总 产 草 量 鲜 重 , 草 层 平均 高 度 是 样 地 
内 多 个 物种 的 高 度 平均 值 。 
2.3 草地 退化 指数 的 遥感 模型 

将 归 一 化 处 理 后 的 植被 覆盖 度 、 草 层 平均 高 
度 、 总 产 草 量 3 个 变量 通过 SPSS 提 取出 1 个 主 成 分 ， 
通过 成 分 矩阵 和 分 指标 的 方差 贡献 率 计 算得 到 综 
合 得 分 模型 中 的 系数 ( 表 5) ,最 后 对 其 进行 归 一 化 
处 理 和 求 平 均值 。 

植被 覆盖 度 . 草 层 平均 高 度 .总 产 草 量 的 权重 
分 别 为 :0.46.0.30 和 0.24( 表 5)。 草 地 退化 指数 
GDI, 计算 公式 如 下 : 

GDI, =0.46 x C,+0.30XH,+0.24xY, (5) 

式 中 : GDI, 表示 地 面 调查 点 位 的 草地 退化 指数 ; C, 
表示 i 群落 标准 化 的 植被 覆盖 度 ; ,表示 i 群落 标 
准 化 的 草 层 平均 高 度 ; Y, 表示 i 群落 标准 化 的 总 产 
草 量 。 
建 模样 地 GDI, 与 6 种 植被 指数 相关 性 最 高 的 
为 NDVILI, 皮 和 尔 逊 相关 性 系数 为 0.743 ,sig<0.01 ,在 
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表 3 20 世 纪 80 年 代 草原 普查 数据 集 
Tab. 3 Grassland census data set in 1980s 
草地 类 草地 型 平均 盖 度 /% “总 产 草 量 /kg.hm3 ”平均 高 度 /em GDI, 
ye fa FE 白 尖 蔓草 ARSE 85 3108.00 14.0 0.52 
鬼 箭 锦 鸡 儿 ZG I ep Be ARS 90 3199.50 25.0 0.78 
线 叶 高 草 90 3237.00 20.0 0.68 
Delta RE AGREE RE , 杂 类 草 90 3031.50 20.0 0.68 
Qe fee ee 、 细 果 蔓 草 PR EE 90 6631.50 17.5 0.75 
Ze My EL 、 杂 类 草 90 7485.00 25.0 0.93 
PRE . 细 果 蔓草 . 杂 类 草 90 3247.50 12.5 0.53 
高 寒 草 原 类 新 疆 曹 草 VA BL 60 1546.50 16.5 0.34 
PAE 65 2173.50 15.0 0.37 
新 好 40 5620.50 12.5 0.26 
新 好 60 4480.50 20.0 0.51 
50 1705.50 10.0 0.15 
KRUSE 85 1770.00 15.0 0.50 
50 1827.00 12.5 0.20 
35 1380.60 20.0 0.23 
35 1380.60 20.0 0.23 
KRUSE AE KEE 85 1914.00 17.5 0.55 
RUZE REFF 85 1914.00 17.5 0.55 
时 黑 麦 . 杂 类 草 85 5319.00 22.5 0.77 
紫花 针 茅 . 冰 草 、 杂 类 草 50 498.00 12.5 0.16 
FRAC Ae ERE HE 50 951.00 17.5 0.27 
紫花 针 茅 .新疆 蔓草 REFF 50 1827.00 12.5 0.20 
KEE REFF kE 50 1827.00 12.5 0.20 
KEE ERE 50 951.00 17.5 0.27 
KEETE OS 45 913.50 17.5 0.24 
紫花 针 茅 、 杂 类 草 50 951.00 17.5 0.27 
座 花 针 茅 、 宽 穗 赖 草 50 1380.60 20.0 0.34 
沼泽 草地 类 HEE RE 100 4128.00 17.5 0.73 
布 顿 大 麦草 、 遍 灯 改 草 100 4323.00 12.5 0.64 


0.01 水 平 上 相关 性 显著 , 优 于 DVI 的 0.694”、RVI 
的 0.568" EVI 的 0.642", SAVI 的 0.637" & TGDVI 
的 -0.435”, 其 中 ,* 代 表 在 0.05 级别 ( 双 尾 ) 相 关 性 
显著 ,** 在 0.01 级 别 ( 双 尾 ) 相 关 性 显著 。 在 GDI, 
为 因 变 量 NDVI 为 自 变 量 的 5 种 回归 模型 (线性 、 指 
数 、 对 数 、 二 次 、 窜 函 数 ) 中 ,线性 函数 模型 显示 出 了 
最 佳 拟 合 ( 尺 =0.553 ) 。 草 地 退化 指数 GDI, 的 计 
算 方程 如 下 所 示 : 
GDI, =0.830xNDVI ,+0.086 
( R° =0.553, RMSE = 0.145 ) 
式 中 : GDI, 表示 草地 退化 指数 的 遥感 拟 合 值 ; 
NDV Lo, 表示 2021 年 的 NDVI 值 。 
为 了 测试 GDI, 回归 模型 遥感 拟 合 值 的 准确 


(6) 


性 ,采用 剩余 30% 的 样 地 点 位 (17 个 样 地 点 位 ) 用 于 
模型 精度 验证 , GDI, 与 GDI, 之 间 RMSE 为 0.140， 
RMSE 值 较 小 ,说 明 反 演 模型 精度 较 高 ,可 用 于 草原 
退化 程度 评价 。 
24 巴 音 布鲁克 草原 的 退化 现状 

相 比 于 20 世 纪 80 年 代 ,2021 年 巴 音 布鲁克 草 
原 退 化 分 级 以 未 退化 为 主 ,未 退化 面积 占 总 面积 的 
比例 为 60.51% ,退化 面积 占 总 面积 的 比例 为 
39.49% , 其 中 重度 退化 面积 .中 度 退 化 面积 和 轻 度 
退化 面积 分 别 占 总 面积 的 比例 为 14.71% 12.55% 和 
12.23%。 

在 3 种 草地 类 型 中 ,高 寒 草 馈 类 草地 面积 最 
大 。 高 寒 草 甸 类 重度 退化 面积 占 该 类 草地 总 面积 
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表 4 2021 年 建 模样 地 数据 集 


Tab.4 Data set of construction points in 2021 


样 地 号 ”平均 盖 度 /% ”总 产 草 量 /(kg:hm”) > 


F 均 高 度 /em GDI, 


001 75 
002 95 
003 100 
004 94 
005 100 
006 80 
007 84 
008 82 
009 60 
010 100 
011 70 
012 70 
013 50 
014 88 
015 80 
016 JI 
017 65 
018 55 
019 65 
020 75 
021 55 
022 55 
023 90 
024 40 
025 85 
026 40 
027 41 
028 40 
029 34 
030 44 
031 92 
032 95 
033 95 
034 84 
035 83 


#5 指标 权重 系数 


3560.00 
4336.00 
6387.00 
2013.00 
10233.00 
2263.33 
2506.00 
2586.00 
1933.00 
1870.00 
1703.00 
2503.00 
2556.00 
3566.00 
1817.00 
2460.00 
1403.33 
1160.00 
776.67 
2070.00 
1210.00 
520.00 
2393.33 
703.33 
1576.67 
710.00 
1980.00 
979.50 
960.00 
864.00 
6853.00 
2660.00 
3240.00 
2482.50 
3720.00 


15.0 
15.0 
20.0 
6.0 
18.0 
7.0 
31.0 
12.0 
12.0 
12.0 
12.0 
22.0 
23.0 
27.0 
16.0 
15.0 
5.7 
5.3 
4.0 
8.7 
5.0 
3.7 
7.7 
4.0 
8.3 
4.3 
6.4 
7.6 
9.6 
4.7 
34.0 
19.0 
9.7 
13.0 
11.5 


Tab. 5 Index weight coefficient 


0.56 
0.70 
0.85 
0.52 
0.84 
0.46 
0.72 
0.53 
0.38 
0.61 
0.43 
0.56 
0.45 
0.74 
0.53 
0.52 
0.33 
0.26 
0.29 
0.44 
0.26 
0.22 
0.53 
0.15 
0.48 
0.15 
0.22 
0.19 
0.18 
0.18 
0.95 
0.67 
0.61 
0.55 
0.57 


植被 覆盖 度 ” 草 层 3 


F 均 高 度 总 产 草 


ta 


20 世 纪 80 年 代 模 型 系数 
2021 年 模型 系数 

20 世 纪 80 年 代 指 标 权 导 
2021 年 指标 权重 
权重 平均 值 


pap 


0.59 
0.59 
0.47 
0.44 
0.46 


0.36 
0.44 
0.29 
0.32 
0.30 


0.30 
0.32 
0.24 
0.24 
0.24 


的 比例 为 21.19% ,大 于 高 寒 草原 类 的 0.67% 和 沼泽 
草地 类 的 14.39% ,高 寒 草 原 类 未 退化 面积 占 该 类 草 
地 总 面积 的 91.54% ,大 于 高 寒 草 向 类 的 49.06% 和 沼 
泽 草地 类 的 41.61%。 高 寒 草 多 类 草地 退化 现状 比 
高 寒 草 原 类 和 沼泽 类 更 为 严重 ,高 寒 草原 类 草地 植 
被 退化 状况 最 轻 。 

由 于 巴 音 布鲁克 草原 山 间 盆 地 地 形 特殊 ,海拔 
相差 大 ,草地 主要 集中 分 布 在 2347~2663 m , 2663~ 
2980 m, 2980~3322 m 的 3 个 海拔 梯度 内 (图 2)。 
1624~2347 m 、3698~4606 m 海 拔 梯度 内 重度 退化 草 
地 面积 占 该 梯度 草地 总 面积 的 比例 分 别 为 63.03%、 
57.32% ,山地 和 平原 地 区 草地 退化 现状 严重 ;2347~ 
2663 m ,2663~2980 m 海 拔 梯度 内 未 退化 草地 面积 
占 该 梯度 草地 总 面积 的 比例 为 75.03% 73.40% ,中 
高 海拔 地 区 尤其 是 盆地 内 草地 退化 状况 较 轻 。 


3.6 

45 | BS 重度 退化 
t 中 度 退化 

2.8 轻 度 退化 


O 未 退化 


2.4 
2.0 
1.6 


面积 /105hm2 


1.2 
0.8 


0.4 


0.0 | 


图 2 不 同 海拔 梯度 退化 面积 
Fig. 2 Degradation area map of different altitude gradients 


2.5 GDI 模 型 的 草原 退化 评价 应 用 
GDI, 模型 可 以 应 用 于 前 期 年 份 的 草原 退化 程 

度 评 价 ,由 于 缺乏 历史 年 份 的 地 面 调查 数据 ,又 不 
能 将 2021 年 GDI, 模型 直接 应 用 于 其 他 年 份 , 故 需 
要 通过 2021 年 NDVI 值 对 2000 年 .2009 年 的 NDVI 
值 做 辐射 配 准 ,提取 2021 年 52 个 样 地 点 位 对 应 的 
2000 年 .2009 年 和 2021 年 的 NDVI 值 ,通过 SPSS 
25.0 完 成 回归 分 析 获 得 辐射 配 准 模型 。 

NDVL =0.901 x NDV Lo — 0.014 

( R =0.682,RMSE = 0.104 ) 


(7) 
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表 6 2000 年 .2009 年 和 2021 年 草原 退化 情况 
Tab. 6 Vegetation degradation in 2000, 2009 and 2021 
1986 一 2000 年 1986 一 2000 年 1986—20094F “1986 一 2009 年 比 1986—2021 年 1986 一 2021 年 
ERSE 面积 /ha 比例 /% 面积 /ha 例 /% 面积 hm 比例 /% 
重度 退化 247967.37 16.47 163884.51 10.89 221458.95 14.71 
中 度 退化 431370.63 28.66 193017.06 12.82 188894.07 12.55 
轻 度 退化 243514.08 16.18 177344.28 11.78 184102.11 12.23 
未 退化 582457.77 38.69 971064.00 64.51 910854.72 60.51 


NDVI; = 0.922 x NDV Loo + 0.070 
( R° =0.680, RMSE =0.104 ) 

根据 公式 (7) 2850 (8 ) XT 2000 4 All 2009 4E ND- 
VI 值 进行 波段 运算 ,将 计算 结果 代入 公式 (6) 计 算 
2 af) GDI, 值 , 最 后 通过 GDI, 的 变化 率 分 级 确定 
2000 年 和 2009 年 草原 退化 程度 ,评价 近 20 a 的 草原 
退化 变化 特征 。 

通过 2000 年 .2009 年 和 2021 年 3 期 退化 结果 明 
显 看 出 ( 表 6),2000 年 以 来 , 巴 音 布鲁克 退化 情况 有 
了 一 定 的 改善 ,2000 一 2009 年 重度 退化 .中 度 退 化 
和 轻 度 退化 面积 占 总 面积 的 比例 分 别 下 降 了 
5.58% 、15.84% .4.4% ,未 退化 面积 占 总 面积 的 比例 
增加 了 25.82% ,退化 情况 得 到 了 明显 的 改善 ; 
2009—2021 年 变化 幅度 不 大 ,重度 退化 、 轻 度 退 化 
面积 占 总 面积 的 比例 提高 了 3.82% .0.45% ,中 度 退 
化 的 面积 占 总 面积 的 比例 下 降 了 0.27% ,未 退化 面 
积 占 总 面积 的 比例 降低 了 4% ,草原 退化 程度 有 轻 
微波 动 ,但 未 退化 面积 仍 大 于 退化 面积 。 


(8) 


3 讨论 


巴 音 布鲁克 草原 海拔 较 高 ,通过 遥感 手段 监测 
时 难以 找到 无 云 干扰 的 多 期 影像 ,长 时 间 序 列 研究 
比较 困难 ,本 文 尽量 保证 了 每 10 a 评 估 1 次 的 逻辑 
性 。 与 其 他 学 者 草地 退化 遥感 评估 的 区 别 在 于 本 
文 形 成 了 20 世 纪 80 年 代 草 原 背 景 值 , 以 GDIL 的 变 
化 率 来 表征 了 退化 程度 。 目 前 存在 3 点 问题 需要 深 
和 人 讨论: 
3.1 草原 退化 遥感 评估 的 误差 分 析 

本 研究 中 GDI, 模型 拟 合 度 为 0.553 ,相关 性 一 
般 , 主 要 有 3 个 方面 的 原因 导致 :(1) 由 于 巴 音 布 鲁 
克 遥 感 影像 云 量 较 大 ,选择 和 地 面 采 样 点 时 间 最 为 
接近 的 Landsat 数据, 二 者 之 间 仍 然 存 在 时 间 异 质 
性 ,时 间 上 不 能 做 到 完全 匹配 ;(2) Landsat 影像 分 


辩 率 为 30 m, 但 样 方 为 1 mx1 m, 空 间 异 质 性 会 导致 
草原 退化 评估 结果 产生 误差 ;(3) 样 方 点 位 过 少 , 过 
拟 合 现象 会 导致 评 佑 结果 产生 误差 。 针 对 草原 退 
化 的 本 底 参 照 系统 , 以 评价 指标 的 最 大 值 “ 或 平均 
值 作为 理想 样 地 参数 也 是 一 种 重要 思路 ,但 会 有 退 
化 评估 程度 过 低 的 问题 ,还 有 学 者 选取 起 始 年 份 作 
为 背景 值 ” ,但 起 始 年 份 草原 的 退化 情况 难以 准确 
把 握 。 为 了 更 接近 草原 未 退化 实际 情况 ,本 研究 的 
巴 音 布鲁克 参照 系统 按照 20 世 纪 80 年 代 草地 类 型 
大 类 赋值 ,将 退化 背景 值 控制 在 合适 范围 内 。 
3.2 GDI 模 型 草原 退化 评估 的 适用 性 

目前 ,应 用 于 草地 生态 系统 的 遥感 应 用 研究 
中 ,存在 定量 反 演 精度 不 高 的 问题 ,本 文 构建 的 
GDI 模型 可 以 解决 草原 退化 定量 评估 的 科学 问题 ; 
并 且 通 过 辐射 配 准 可 以 完成 GDI 模 型 在 前 期 年 份 
的 应 用 ,填补 由 于 历史 数据 缺失 导致 连续 时 间 序 列 
研究 的 空白 ,时 间 尺 度 上 模型 更 具 优 势 , 可 以 用 于 
指导 草原 生态 系统 保护 的 数据 基础 ,对 于 长 时 间 序 
列 草 原 退化 评估 该 模型 普 适 性 更 高 。 并 且 通 过 引 
和 人 权重 可 以 实现 多 指标 的 综合 ,提高 遥感 草原 退化 
评估 的 效率 ,满足 草原 生态 系统 管理 的 应 用 需要 。 

2000 年 未 退化 面积 占 总 面积 的 比例 为 38.69%， 
2009 年 未 退化 面积 占 总 面积 的 比例 为 64.51% ,2021 
年 未 退化 面积 占 总 面积 的 比例 为 60.51% ,本文 结 
与 焦 阿 永 等 3 通过 NDVI 值 计算 得 出 的 2000 一 2019 
年 北 疆 地 区 草地 禾 盖 显著 增加 的 结论 相符 ;但 其 中 
2000 一 2009 年 草原 退化 程度 相差 较 大 ,新 疆 草地 植 
被 - 净 初 级 生产 力 与 降水 呈正 相关 关系 ,通过 辐 
射 配 准 方法 完成 前 期 年 份 的 草原 退化 评估 存在 气 
象 背 景 值 不 统一 的 误差 ,优先 完成 研究 区 不 同 降水 
年 份 的 分 类 ,可 以 提高 模型 的 精度 ,后 期 工作 需 
要 把 巴 音 布鲁克 草原 分 为 降水 平年 FARF, 
比较 相同 降水 背景 下 的 草原 变化 情况 。 


cf 


202305.00262v1 


chinaXiv 


3.3 草原 退化 的 影响 因素 

植被 长 势 好 坏 与 气候 因素 密切 相关 ,Fu 等 * 研 
究 显示 ,1961 一 2009 年 新 疆 天 山 开 都 河流 域 气温 以 
0.167 'C…(10a) 的 趋势 增加 ,降水 以 0.167 °C -(10a)' 
的 趋势 增加 ; 刘 漂 漂 等 中 对 巴 音 布 鲁 克 气候 的 研究 
结果 表明 ,20 世纪 90 年 代 中 期 以 后 , 巴 音 布鲁克 气 
温 升 高 .降水 增多 ,因此 本 文 研 究 时 间 段 内 的 
1986 一 2015 年 巴 音 布鲁克 草原 气候 总 体 呈 现 “ 暧 湿 
化 ”; Gang 等 “人 研究 表明 ,新疆 1981 一 2010 年 的 草 
地 NPP 同 降水 和 气温 呈正 相关 关系 , 巴 音 布鲁克 草 
原 温度 回暖 和 降水 增加 趋势 使 得 当地 草地 植被 得 
以 恢复 ,这 与 赵 鹏 等 汪 研 究 发 现 ,1982 一 2015 FA 
候 要 素 导 致 新 疆 草地 NPP 增 加 的 结果 相似 ,与 本 文 
研究 得 出 2000 一 2009 年 草原 退化 程度 明显 改善 的 
结果 相符 ,而 对 于 2009 一 2021 年 草原 退化 程度 轻微 
波动 的 结果 , 姚 俊 强 等 研究 结果 指出 ,20 世 纪 90 
年 代 之 后 新 疆 气 候 多 呈 暧 湿 配 置 ,但 在 2022 年 时 段 
内 有 所 减缓 ， 暖 湿 化 "趋势 减缓 导致 了 2021 年 草原 
植被 状态 的 非 线性 增加 结果 。 此 外 ,人 类 活动 对 草 
原 的 影响 复杂 多 变 , 曹 永 香 等 ”通过 主 成 分 分 析 法 
研究 发 现 , 人 类 活动 是 策 勒 绿洲 -沙漠 过 渡 带 植被 
变化 的 主导 影响 因素 。 除 了 牲畜 量 的 增加 导致 草 
原 退 化 外 , 退 牧 还 草 工程 下 的 围栏 封 育 、 季 节 性 休 
牧 、 划 区 轮 牧 、 牧 革 补 播 等 是 目前 人 类 保护 草原 
生态 的 主要 措施 ,人 类 活动 在 草原 退化 程度 改善 中 
起 到 了 多 大 的 作用 需要 选取 合适 的 科学 方法 进 一 
步 分 析 。 


4 结 论 


本 文通 过 构建 GDI. 模型 和 Landsat 影 像 对 巴 音 
布鲁克 草原 多 年 退化 程度 进行 评价 ,提高 了 草原 退 
化 遥感 监测 效率 , 主要 结论 有 : 

(1) GDI, 综合 指标 为 植被 覆盖 度 、 草 层 平均 高 
度 、 总 产 草 量 的 加 权 平 均 , 引 入 该 指标 可 以 提高 草 
原 退 化 遥感 评价 的 效率 ; GDI, 与 NDVI 的 相关 性 最 
好 , 且 GDI, 模型 方程 的 拟 合 度 为 0.553, RMSE 为 
0.140 ,适用 于 巴 音 布鲁克 草原 退化 程度 遥感 评估 。 

(2) 2021 年 巴 音 布鲁克 草原 未 退化 面积 占 总 面 
积 的 比例 为 60.51% ,高 寒 草 多 类 草地 退化 最 严重 
高 寒 草 原 类 草地 退化 状况 最 轻 , 空 间 分 布 上 表现 为 
由 盆地 向 山地 退化 加 重 的 趋势 。 
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(3) 利用 辐射 配 准 方法 对 前 期 年 份 草原 退化 程 
度 进行 评估 ,研究 的 时 间 尺 度 更 长 ,结果 表明 
2000 一 2009 年 巴 音 布鲁克 草原 退化 程度 明显 改善 ， 
2009 一 2021 年 草原 退化 程度 轻微 波动 。 
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Evaluation of the degree of degradation of Xinjiang Tianshan 
Bayinbuluk grassland in 35 years 
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Abstract: Many studies on grassland ecosystem degradation have been performed using vegetation coverage or 
productivity indexes for evaluation. However, it is difficult to comprehensively evaluate different degrees of 
grassland degradation using a single evaluation index. Taking Bayinbuluk grassland in the Xinjiang Tianshan 
Mountains as a research object, a remote sensing method for evaluating grassland degradation based on 
standardized processing sub- index coupling was proposed. The grassland vegetation coverage, average grass 
layer height, and total grass yield were selected to determine the weight of the three indexes by principal 
component analysis. The Min- Max standardized method was introduced to construct the grassland degradation 
index (GDI). Finally, the degree of degradation of Bayinbuluk grassland in the Xinjiang Tianshan Mountains was 
determined through Landsat image inversion and reasonable classification of the rate of change of grassland 
degradation index from 1986 to 2021. The results showed that GDI, has the best correlation with NDVI. In 2021, 
the proportion of undegraded area of Bayinbuluk grassland relative to the total area was 60.51%. The degree of 
degradation of different grassland types showed clear differences. The spatial distribution showed a trend of basin 
to mountain degradation. The GDI model could be applied to other years through the radiation registration 
method. The degree of degradation of Bayinbuluk grassland significantly improved from 2000 to 2009 and 
slightly fluctuated from 2009 to 2021. The results of this research provide robust data support and a theoretical 
basis for guiding evaluation of the degree of degradation of Bayinbuluk grassland and protecting the grassland 
ecosystem. 

Keywords: grassland degradation; Landsat image; Min-Max standardization; grassland degradation index; 
regression model; Bayinbuluk grassland 


